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We show that several classes of mixed quantum states in finite-dimensional Hilbert spaces which 
can be characterized as being, in some respect, 'most classical' can be described and analyzed in a 
unified way. Among the states we consider are separable states of distinguishable particles, uncor- 
related states of indistinguishable fermions and bosons, as well as mixed spin states decomposable 
into probabilistic mixtures of pure coherent states. The latter were the subject of the recent paper 
by Giraud et. al. [Jl, who showed that in the lowest-dimensional, nontrivial case of spin 1, each such 
state can be decomposed into a mixture of eight pure states. Using our method we prove that in 
fact four pure states always suffice. 

PACS numbers: 03.65.-w, 03.67.Mn, 03.65.Fd 



The notion of being the 'most classical' among all quantum states of a given system is admittedly vague, and strongly 
dependent on the particular quantum features we would like to approximate on the classical level - 'nonclassical' 
properties of quantum systems usually encompass a variety of phenomena. 

In the situation that we will be concerned with here, there is a Lie group of preferred observables acting on the 
system. It may be the group SU{N) x SU{M) of local unitaries acting on a composite system, or the SU{2) group 
one can conveniently realize using linear optics. Then the group action divides the set of states into disjoint orbits, 
and there will be a special orbit which is in itself a symplectic and indeed a Kahler manifold. It can serve as a classical 
phase space, and the states in this orbit are our 'most classical' states. They are separable states in the first example, 
and coherent spin states in the second. In general they are the generalized coherent states defined as the orbit 
in the representation space of the highest weight vector for a particular irreducible representation of the symmetry 
group of the system. Generalized coherent states also minimize an uncertainty relation 3], or more precisely they 
minimize the total variance of a state ip, 



where the {Xi} form an orthonormal basis in the Lie algebra of the relevant group. Note that this is independent of 
the choice of the basis {-'^i}. It is pleasing that both definitions of 'most classical' coincide, and moreover, that they 
coincide with the view that separable states of a composite system are the most classical ones, in the sense that they 
lack the quantum correlations among subsystems responsible for any EPR-like phenomena. 

In the following we call (for shortness) those states 'classical' which are 'most classical' in the above senses. At the 
other end of the spectrum we can identify the states in that special orbit on which the symplectic form vanishes as 
being the 'most non-classical' states. The reason why there is such an orbit can be seen by complexifying the group; 
for composite system this coincides with the maximally entangled states as conventionally understood |J] . 

Investigations concerning role of pure (generalized) coherent states in exhibiting links between quantum and classical 
level have a long history, especially for such problems as quantization and classical limits (see eg. Q)- The connections 
between coherent states and entanglement theory was recently pointed by Klyachko 0,0 (see also 0, Bll])- 

The notion of classicality can be extended to mixed states. On a formal level we do it by defining 'classical' mixed 
states as those which can be decomposed in a probabilistic mixture (i.e. a convex linear combination) of pure classical 
states (here and in the following we identify pure states, i.e. vectors in a Hilbert space, or rather its projectivisation, 
with one-dimensional projection on their directions). For composite systems it leads to the popular definition of 
separable (or equivalently non-entangled) states, which can be also characterized on the operational level as those 
which can be obtained from an initially uncorrelated state of the whole system by local operations, i.e. quantum 
operations performed only on subsystems, accompanied by possible exchange of classical information between them 
(e.g. about results of local measurement on which further local operations can be conditioned) [To| . The general 
practical criteria useful for unambiguous discrimination of mixed separable states among all states of a composite 
quantum system are not known, but in low-dimensional cases there are several methods (e.g. based on determining 
the so called concurrence of a state) to achieve the goal. 



I. INTRODUCTION 




(1) 
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Equating absence of genuine quantum correlations with the property of being a product state (simple tensor) lacks 
sense in the case of indistinguishable particles. Indistinguishability forces (anti)symmetrization of wavefunctions, in 
other words states are vectors not in tensor products of Hilbert spaces but in their (anti)symmetrizations. With 
the exception bosons occupying the same state, no wave function has the form of a simple tensor. It is reasonable 
to call 'nonclassicaP quantum correlations which go beyond those stemming solely from (anti)symmetrization, and 
consequently states which lack such correlations are treated as quantum uncorrelated or the best candidates for being 
'most classical'. Methods of investigating and characterizing such states were developed in [ll| and [l^ . 

For spin-coherent states, when the group is SU{2), the classical mixed states — as is customary in quantum optics — 
are characterized as having a positive P-representation. The definition is equivalent to decomposability of a state into 
a probabilistic mixture of pure coherent states. For the lowest-dimensional nontrivial case (spin 1) Giraud et. al. [l| 
give a necessary and sufficient condition for classicality of the state. Their proof of the correctness of the criterion is 
constructive — it gives explicitly a decomposition into eight coherent states. 

The observation that all above mentioned cases are instances of the general construction of coherent states allows 
to analyze them in a unified manner. One may thus use methods known from the theory of coherent states to 
the separability problem, or vice versa, apply techniques from entanglement theory to exhibit some new features of 
generalized coherent states. We will mainly follow the latter path — a concrete aim is to refine one of the results 
of [l| concerning the cardinality of the optimal decomposition of a mixed spin-coherent state into a mixture of pure 
coherent states. Along the way we present a unified treatment of correlations in two-component quantum systems 
in terms of bilinear observables. We start by recalling the characterization of entanglement by the (pre)concurrence, 
defined originally by Wootters [ll] and elaborated further by Uhlmann In this approach one quantifies the 

entanglement of a (pure) state by the expectation value of some antilinear operator and extends the measure to 
mixed states by the 'convex roof construction, i.e. by decomposing a mixed state into a probabilistic mixture of pure 
states and minimizing over all possible decompositions the sum of concurrences of the components. As observed by 
Uhlmann antilinear operators 'are intrinsically nonlocal', and as such are good candidates for a tool probing 
nonlocal properties like entanglement. On the other hand, they do not correspond to physical observables, hence their 
expectation values can not be directly measured. Elsewhere p^ . [l6j we argued that it might be more convenient to 
use bilinear Hermitian operators to quantify entanglement — such an approach is efficient in higher dimensions and 
provides measures of entanglement accessible directly in experiments [17| . Remarkably the same approach enables 
also a unified treatment of the different facets of 'classicality' enumerated above. 



We start with an elucidation of links between generalized concurrences defined by Uhlmann and bilinear operators. 
As already advertised in the previous section, both provide tools for characterizing separability and entanglement as 
we will show in Section HVl 

In the following TC will denote a finite-dimensional complex Hilbert space and is an orthogonal basis in 

Ti,. Let ip Ti and Cq{iP) — the generalized concurrence in the sense of Uhlmann [14| — be defined via an antiunitary 
operator 9, i.e. 



Each antiunitary O can be written as O = Tq K, where Tq is unitary and K is the complex conjugation in some fixed 
basis (say the chosen one {\ei)}iLi)^ i-e. 



II. GENERALIZED CONCURRENCES 



Ce(^) = |(^|e|^|)|. 



(2) 



Ce{i^) = mTe\r)\ 



(3) 



Let us define A e End{H (K" H), 



A= (I® A)|$)($|, 



(4) 



where |<I>) is the maximally entangled state in TC (E)'H, 



N 




(5) 



and for an arbitrary p e End{H) 



A(p) = TepT, 



0' 



(6) 
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i.e. A is the image under the Jamiolkowski isomorphism of A which, in turn, is a completely-positive unitary map on 
End{H (8) H). Since A is completely-positive, A is non- negatively definite [Tol. [isl [Tot. 

The announced connection between the bilinear operator A and the generalized concurrence Cq is given by the 
following 

Lemma 1 

CeW = (7) 

Proof 

From the definition of |$) and A 

A= (l®A)^|e,)(e,|®|e,)(e,| =^|e,)(e,|^Te|e,)(e,|4. (8) 

Hence, for = J2k'^k\ek), 

(V- V'l^lV' ® V') = ^ clciCmCn{ek\et){ej\em}{ei\T0\ei){ej\Tl\e,i) = ^ 4q*c™c„4i^jm(ei|re|ei)(ej|r^|e„) 

ijklmn ijklrnn 

= 5]c*crc,c„(ez|re|e,)(e,|4|e„> === {i^\Te\r) {r K\^) = WTe\r){i'\Te\rr = MTe\r)? = Ceii^f ■ (9) 

ijln 

□ 

Remark 1 We have 

TiiA ^ I, (10) 

and 

Tt2A = I, (11) 

where Tri^2 denotes the partial trace over the first or the second copy ofTi.. Indeed, from the unitarity ofTo, 

TriA = ^(efc|e,)(e,|efc)re|e,)(e,|4 = ^ 4z'5,fcre|e,)(e, |4 = re(5] |e,)(e,|)4 = Tq^ = I, (12) 

-ij ijk i 

and 



Tr^A = |e»>(e,l(efc|Te|e.)(e,|4|efc) = ^ |e,)(e,|(4)^.,(Te),^ = ^ |e.)(e,|(4re 

ijfe ijk ij 

= ^|e»)(e,|<5,,=^|e,)(e,| =L (13) 

ij i 

Observe now that we may proceed in the opposite way, i.e. define the generalized concurrence as in ([7]), 

CaW = (V- ® ® V')'/', (14) 

where A is some non- negatively defined linear operator oii Ti ® Ti with the property (|10p . Then we define A G 
End{End{'H)) as the image of A under the inverse of the Jamiolkowski isomorphism, i.e., 

A(p) =tri((p*®I)A), (15) 

for p £ End{H). Observe that from ([10|) and (fT5| . 

A(I) =tri(A) =1 (16) 
Since A is non- negatively defined, A is completely positive, hence it has a Kraus decomposition [lOllTgj. 

s 

A{p)^J2^^P^^. (17) 
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If A is such that s = 1, then Ap = TipT{ and from (HH) TiT/ = I, i.e. Ti is unitary. Then we have 

CA{^P)^MTi\r)\^CeW, (18) 

where Q = TiK. To prove the last statement it is enough to check that (fT5|) gives indeed the inverse of the Jamiolkowski 
isomorphism ([¥]) and perform the calculations in the proof of Lemma [T] in the reverse direction. 

It seems to be more appropriate to take (|14l) with the condition pO)) as the definition of the generalized concurrence 
(see Generically the non-negative definite matrix A has the spectral decomposition 

s s 

A ='^iya\Va){Va\ =^^\'Wa){Wa\, (19) 
a—1 a—1 

with more then one nonvanishing eigenvalues Va- Here \wa) — y/ly^\va) are the subnormalized eigenvectors of A. It 
means that A defined by ([15]) is given by pT)) with 

T„= (20) 

Only exceptionally (eg. in low-dimensional cases) s = 1 and we can write Ca in terms of an antiunitary operator O. 
Performing the same calculations as in ^ for the general case ([TT]) we obtain: 

s s 

(V'® vi^iV'®^) = 5](V'ir„i^*)(viTiiv) = MTa\r)\\ (21) 

a—1 Q— 1 

and by polarization, 

s 

(V'2®V'4iAiV'i®V'3) = Y.^HTam){ri\Tim- (22) 



III. MIXED STATES 

The concept of concurrence can be extended to mixed states. The idea is based on the following observation [23| 
(see also Bill). 

Let E be the set of all extreme points of a compact convex set in a finite dimensional real vector space V. For 
every non-negative function f : E ^ we may define its extension fx ■ K — > K_|_ by 

fK{x)= inf ^p./(x,) (23) 

where the infimum is taken with respect to all expressions of x in the form of convex combinations of points Xi from 
E. Let now Eq be a compact subset of E with the convex hull Ko — conv{E{j) G K. If / is continuous and vanishes 
exactly on Eq, then the function fx is convex on K and vanishes exactly on Kq. In our cases K is the set of all 
states, E — the set of pure states, and Eq — the set of pure 'classical' states. Observe that due to homogeneity of 
the generalized concurrence, C^(q;'0) = |apCA(V')i we can consider only the decompositions into sums of rank-one 
operators, defining thus for an arbitrary state p, 

K 



CA{p)^mm\CA{<Pk), (24) 



where the minimum is taken over all decompositions of p into a sum of rank-one operators, 

K 



fc=i 



We took advantage of the finite dimensionality of the Hilbert space and substituted the minimum for the the infimum 
operation [2^ — in the finite-dimensional case there exists always an optimal decomposition, i.e. one which minimizes 

CAip). 
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A particular example of ([25)1 can be obtained from the spectral decomposition of p, 

r 

P^^P'Mi){ri,\, p\m) ^ Pilvi) , {Vi\Vj)^Sij, r^rankp, (26) 

i=l 

by subnormalizing the eigenvectors, 

r 

\C.) = V^^\V^), P^T.\^^m\■ (27) 

i=l 

Any other decomposition (PS)) can be obtained from (^7)1 with the help of a partial isometry [2^ . 

r 

\<Pk)^Y.^kM^), k=l,...,K; VW^I. (28) 

Hence, 

Ca{p) = min J2 C A{(t>k) ='aiin^{(f>k<»(l)k\A\(pk^<pky^^ 

k k 

= min^ ® ® 0™))'^' 

1/2 



where Tq, are r x r matrices, 

(Ta)., = (e.lT^alf*), (30) 

and the minimum is taken over all partial isometrics V . This can be a starting point for various estimations of Ca{p) 
in the spirit of [l^ and p^ . 

In the case when s = 1 (so the generalized concurrence for pure states can be expressed in terms of an antiunitary 
operator (O), the expression ([29]) reduces to 

Ca{p) = min^ \{V*TV^)kkl n, = (^l^iie*). (31) 

A; 

If Ti is symmetric, Ti — , (as it happens in all the cases we consider - see below), then = TiK is an antiunitary 
conjugation, i.e. 9^ = 1, and Ca{p) can be explicitly calculated, 

r 

Ca(p) = max |o,^i - ^/ijj, (32) 

where /i^ are singular values of t in decreasing order. A short calculation shows that fii are equal to the square roots 
of the eigenvalues of pp, where p = QpQ. The eigenvalues of i? := pp are positive since R is similar to a Hermitian 
positive-definite matrix R = p^^^R'p^^^^, 

R' = p^/^QpQp^'^ = p'I^T^p*TIp^I^ = (pi/2Tip*i/2)(^*i/2y*^i/2) ^ (^i/2y^^*i/2)(^i/2j.^^,i/2)t^ (33) 

since p = p^ and Tf = Ti i.e. = Tj*. 

The proof of ([5^ was given by Uhlmann (see also [l^) who also showed that the optimal decomposition may 
be constructed out of 2"+^ vectors, where 2" < TV < 2"+^. This is a generalization of the Wootters construction for 
the concurrence of two qubit states [l^ . 

In the following we will call a (pure or mixed) state 'classical' if an appropriate generalized concurrence (j24p vanishes. 

From the arguments in the last two sections it should be clear that although vanishing of an appropriate Ca{p) 
can characterize 'classical' (eg. separable) states in an arbitrary finite dimension, it is only rarely that its calculation 
can be reduced to (|32p via some antilinear Q. Such situation happens only in low-dimensional cases and this is the 
reason why e.g. Wootters' construction gives explicit results only in the two-qubit case. In the next section we 
will give more examples supporting this observation. 
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IV. EXAMPLES: SEPARABLE AND UNCORRELATED STATES OF COMPOSITE SYSTEMS 
A. Entanglement in two-partite systems of distinguishable particles 



The Hilbert space TC is the tensor product of the Hilbert spaces of the subsystems, H = Hi (8) 7^2 of dimensions 
A''i and N2. In this case A is proportional to the projection on Tiu ■— {Tii A Hi) {H2 A H2) This is a direct 
consequence of the fact that for a pure state \iJj) the quantity 



(34) 



where Pr is the reduced (to one of the subsystems) density matrix, vanishes only for separable {ip), since it is only 
in this case that the reduced density matrix is pure and the trace of its square equals one. It is now a matter of 
straightforward calculations to establish that (l34l) is proportional to (-0 (8) V'l^lV' ® ip) where A is the projection on 

'^12- 

Let us briefly rederive this result using Schmidt decomposition, since later we would like to stress analogies with 
other examples of 'classical' systems. An arbitrary two-partite pure state can be written as 



(35) 



where {|ei)}|iJi and are some orthonormal bases in Hi and H2 and c some complex matrix. An arbitrary 

complex matrix c can be transformed to a diagonal one with non-negative entries by multiplication from left and 
right by unitary matrices, c 1-^ UcV, which amounts to local unitary changes of bases, \ei) = 12k^kiWk) : 

\fj) = yji\fi)- Hence, upon an appropriate choice of U and V we obtain tp in its Schmidt form 



k 



(36) 



where the last condition stems from the normalization Tr|'0)(-!/)| = 1. The state \ip) is a simple tensor (i.e. is separable) 
if and only if only one among Ik docs not vanish. It is now easy to show that this condition is equivalent to the 
vanishing of C given by ([34]) . 

The classical pure states are thus the separable ones. Hence the classical mixed states are those which are decom- 
posable into pure separable states, i.e. (mixed) separable states. 

For two qubits, A^i ~ 2 — N2, the matrix A reads 



A 











































1 



-10 





-10 





1 












-10 





1 




1 



-10 















-10 





1 




1 



-10 















10 





-10 





-10 





10 









(37) 
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It has only one nonvanishing eigenvalue (equal to 4) with the eigenvector 

\wo) = [000 -1 00100100 -1 000]'^. 
The corresponding matrix T reads 



(38) 



T = 



■ 








-1 - 








1 








1 








. -1 








0. 



(39) 



and the concurrence is given by the Wootters formula, 

4 

C{p) = max |o, Ail - X!^^ }' 



(40) 



3=2 



Here fij are the singular values of Tij — {^i\T\£^*) and i = 1, . . . ,4 — the subnormalized eigenvectors of p, or 
equivalently, pi are the square roots of the eigenvalues of pp, where p — Tp*T . The optimal decomposition can be 
constructed out of four vectors, which is also clear from Wootters' original construction. 



B. Entanglement in two-fermion systems 



Correlations in systems of two fermions were investigated in [ll|, \\^ , where all relevant definitions and proofs can 
be found. Here we only briefly review the most important findings. In this case Ti is the antisymmetric part of the 
tensor product of two copies of the single-particle Hilbert space TL2K of an even dimension 2K = 25* + 1, where S is 
the spin of each particle, Ti = A{'H2K ^'H2k), hence H has dimension N = 2K{2K — l)/2 = S{2S + 1). An arbitrary 
pure state can be represented in the form 



2S+1 



1^) = E ^^.-///jio)' 



(41) 



where are fermionic creation operators and |0) — the vacuum state and w — a complex antisymmetric ma- 
trix. A unitary transformation U of the single particle space TCs leads to the transformation fj i—> ^^-^ Ujifj and, 
consequently. 



w I-* UwU 



(42) 



The analogue of the Schmidt decomposition is now provided by the theorem stating that an arbitrary, complex matrix 
w can be brought by an appropriate transformation (|42p to a block-diagonal, canonical form. 



w 



diag[Zi, . . . , Zr, Zo] , Zi 



z, 

-z,, 



(43) 



with Zi ^ and Zq - a. null matrix [23|. The number r of non- vanishing 2x2 blocks is called the Slater rank of lijj) 

- . . . . 

Pure states with the minimal (i.e. equal to one) Slater rank exhibit the minimal allowable quantum correlations 
so they are the candidates for the 'classical' states. Consequently, mixed states are 'classical' when they can 
be decomposed into a convex combination of pure states with Slater rank one. The maximal Slater rank of the 
components of a pure state decomposition of a mixed state p minimized over all possible decomposition is called the 
Slater number, by analogy to the Schmidt number defined in a similar manner for distinguishable particles. Using 
this notion we identify 'classical' mixed states with those of the Slater number one. 

In order to characterize the 'classical' states we make use of a lemma proved in p^ . It provides a general criterion 
for a pure state in the form (|4ip to have a Slater rank not exceeding a prescribed value. In particular it states that a 
two fermion state inTl — AiTi^K ® 'H2K) has the Slater rank one if and only if for all 1 < ai < • • • < a2{K-2) < 



2K 



(44) 
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where e'i - *2it jg the totally antisymmetric unit tensor in Ti.2K, which is obviously equivalent to 

2 



E 



l<ai < ■ ■ ■<» 



2{K -2)<2K 



2K 

wijWkis 

i.j,k,l — l 



ijklax...a2{K-2) 



= 0. 



(45) 



Each of the terms under the absolute value sign is the PfafRan [30| of the 4x4 matrix obtained from w by deleting 
the rows an columns with numbers ai, a2, . . . , a2(x-2)- It is a matter of simple calculations that the Pfaflian Pf(X) 
of a 4 X 4 matrix X reads 



and 



where |x) is the six-dimensional vector 



|pf(x)r = (x«)x|^p/|x(g)x), 



^12 

X34 



(46) 
(47) 



and, in the standard basis |ei) 



lee; m 



(48) 



A 



Pf 



with 



T 



E 










)(e 














" 














1 " 














-1 














1 














1 














-1 














1 


















(49) 



(50) 



Substituting for each term in the sum in (I45[) an appropriate expression (j47p , involving in each summand a different 
set of entries of w with the corresponding matrix Apf, we rewrite the condition (|45p 

C^(?/)) (w(8) w|A|w(8) v^f) = 0. (51) 

Here A is a 5^(25* + 1)^ x 5^(25 + 1)^ matrix acting in 7i 7i, built up from the submatrices Apf of each summand, 
and |w) is the 25+1 dimensional vector of independent entries of the matrix w representing in this way the state 
vector \ip) (I4ip . The condition (|5ip constitutes the appropriate non-entanglement criterion for pure fermionic states 
in the form of a bilinear generalized concurrence (|14p . 

For the lowest-dimensional, nontrivial case, i.e. S = 3/2, yl is a 36 x 36 matrix (skipped to save the place) with 
only one non-vanishing eigenvalue. The corresponding eigenvector has only six non-vanishing elements and the 
corresponding matrix T ([20|) is given by (I50p . as it is also clear from ^ and (|49|). The result is equivalent to the one 



of [141 modulo a change of basis in C®: \e[) — |ei), jej) — |e2), leg) = (jea) -|- |e4))/\/2, \e'^) ~ jes), jcg) — leg), and 
|e^) = (les) - |e4))/V2. ^ ^ ^ . . 

The generalized concurrence Ca{p) for mixed states is the Slater correlation measure |ll| and is given by 



6 

Csi(p) = max|o,/^i - } 



J=2 



(52) 



The definition of /ii is analogous to the one used in the preceding example, the only changes are in the dimensionality 
and the definition of T. As previously fii are equal to the square roots of the eigenvalues of pp, with p = Tp*T, which 
coincides with the results of [ll|, [l^ . From the result of Uhlmann we infer that the optimal decomposition can be 
achieved with 8 vectors. 
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C. Entanglement in two-boson systems 



Quantum correlations in bosonic systems were thoroughly investigated in 12]. The line of thought follows mutatis 
mutandis considerations in the fermionic case. For two particles the Hilbert space is the symmetric part of the tensor 
product of the single-particle space, H = S{Hm ^ Hm), dimTlM = M, dimH = M{M + l)/2. A pure state can be 
written in the form 



M 



(53) 



with bosonic creation operators 5| and a symmetric complex matrix v transforming upon unitary map U in the single 
particle state according to 



'YliiUijhj and, consequently. 



As previously, an appropriate theorem from the linear algebra p4 
to its diagonal, canonical form 



(54) 

provides a possibility of using (|54p to transform v 



diag(zi, 



,0,...,0), z,;^0. 



(55) 



The number r of nonvanishing Zi is dubbed bosonic Slater rank, and states with the minimal r = 1 exhibit the 
minimal possible amount of purely quantum correlations [l^ . Accordingly these are our 'classical' pure states, 
whereas 'classical' mixed states are those which can be decomposed into a convex combination of pure states with the 
Slater rank equal to one. As in the previously considered cases one can define the bosonic Slater number of a mixed 
state as a minimum over all convex decompositions into pure states of the maximal Slater rank among the members 
of a decomposition. In this terminology, a mixed is 'classical' if and only if its bosonic Slater number equals one. 

As in the case of fermions there exists a bilinear characterization of pure states with the bosonic Slater rank equal 
to one [131 . For them 



M 

ijkl—1 



ikai...aM-2 ^jlai...aM-2 



(56) 



for all 1 < ai < • • • < aM-2 < M, which can be transformed to the desired form (fT4|) along, essentially, the same 
lines as in the fermionic case. To this end we rewrite ([55)) in the form 



E 

l<ai<---<Qjj_2<A^ 



M 

E ^ij'^klS 



ika\ . . .aA 



= 0. 



(57) 



To cut short the connection to the conventions used in [12!| it is convenient to represent the vector 1-0) (j53p in 



Vl2 



(58) 



of the independent entries of the symmetric matrix v, i.e. Vij, i < j, scaling, however, the diagonal entries by l/%/2, 
v'jj = Vjj/V2. The way of representing 1-0) in terms of v is, obviously, a matter of convenience, as long as linearity 
of the representation is observed. The chosen one, besides being in accordance with fT5|, has the advantage of giving 
the same 'weight' to diagonal and off-diagonal entries of v. 

Each summand under the absolute value sign in (j57p involves only a 2 x 2 submatrix of u. 



v{a, b) 



Vaa Vab 
Vab Vbb 



l<a,b<M, 



(59) 



and equals, up to an irrelevant sign, 2va 
(v ® v|Af,os|v ® v), where 



iVbb 



-"aa^bb 



This, in turn, can be rewritten in the form 
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|v) = |v(a,b)) = 



Vab 
V'bb 



A, 



bos 



\e,){e,\(g,T\e,){ej\T\ T = 



and {|ei), |e2), jes)} is the standard basis in 
rewritten as 



1- 

0-10 , (60) 
10 0. 

Again, the non-entanglement condition for pure states ([57]) can be 



3 

E 







(8) V 


\A\v (8) 


v) =0, 


by combining the appropriate 6x6 matrices At,os into a (M(M + l)/2)^ 
In the lowest-dimensional, nontrivial case M = 2, iV = 4 we find [3l| 
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(61) 



(62) 



The matrix A has only one non-vanishing eigenvalue with the eigenvector 

[ 1 -1 1 ]"^, 



(63) 



to which there corresponds via (j20p the matrix T given in (j60p . In a perfect analogy with the previous examples the 
generalized concurrence Ca{p) is given in terms of the square roots of the eigenvalues of pTp*T, 



Ca{p) = max|o,/ii - ^Mj}; 

J=2 



(64) 



and coincides with the hosonic correlation measure jl2l |. The optimal decomposition involves at most 4 vectors. 



V. SPIN COHERENT STATES 



Classical spin states were defined in [l| as those which can be decomposed into a probabilistic mixture of pure 
spin coherent states. In order to mimic the previous construction we will use the following characterization of pure 
coherent states, following directly from their minimal-uncertainty property Q. 

Theorem 1 A spin-S state {"tp) £ Ti. is coherent if and only if 

{Li ® Li + L2® L2 + Lz® £3)!^' 8) -0) = 5'^ 1-0 ® -0), (65) 

where Li are operators of the spin components in the spin S representation of the rotation group. 

This is obviously equivalent to 

{tp®'il:\A\ij®'il;) A = I - Lx® Li- L2® L2~ L:^® L^. (66) 

For the lowest-dimensional nontrivial case S* = 1 we choose 



(67) 





- 





0" 




- 





i ' 




■ 


—i 


0" 


Li = 








—i 


, L2 = 











, ^3 = 


i 










. 


i 


0. 




. —i 





0. 
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Then 



A = 



1 

1 










1 





1 

1 




1 

1 



1 
1 













1 



1 









1 



1 




1 
-1 



1 
1 



(68) 



The matrix A has eigenvalues 3, 2, 2, 2, 0, 0, 0, 0, 0, hence, at first sight, the previous construction fails. But by inspec- 
tion of the eigenvectors we find that the corresponding matrices Ta read 



(69) 



hence, due to the antisymmetry of Ti, i = 2,3, 4, we have (V'lTil^*) ~ for i ^ 2,3, 4, and an arbitrary 3-dim vector 
IV'). Thus (cf. (Ell)) 



and 



CaIp) 



-3 

:{0,/ii-Xl^j}' 



(70) 



(71) 



J=2 



where fii are the singular values of = and |^i) are the subnormalized eigenvectors of p (or, equivalently, fii 

are the roots of the eigenvalues of i? pp*)- The optimal decomposition can be achieved with at most four vectors. 

The criterion of classicality based on (j7ip can be formulated directly in terms of orthogonal invariants (traces of 
powers) of R. To this end consider a 3 x 3 density matrix i.e. p = p\ Trp — 1 and p — non- negatively definite. In 
the following we assume that p is nondegenerate, det p > 0, we also will suplement this by some other (see below) 
non-degeneracy conditions (all degenerate cases, in general simpler in treatment, can be considered along the same 
lines). 

We decompose p into the real and imaginary part, p — pb. + ipi- The hermiticity of p implies — pn, pj — —pi- 
The real part p^ can be thus diagonalized by an orthogonal transformation which leaves the antisymmetry of pi 
unaltered. After that p takes the form 



P = 



W3 A2 ~ivi 
-iv2 ivi A3 



(72) 



Since p is non-negatively definite we have A^ > 0. In the following we assume a further non-degeneracy condition 
Aj > 0. We have also Ai -t- A2 + A3 = Trp = 1. 
The matrix R reads now: 



R = pp* 



V2 



V3 



V1V2 + ivsiXi ~ X2) V1V3 + iv2{X3 - Xi) 



U1U3 + W2(A3 - Ai) A2 



V3 



V2V3 + ivi{X2 - A3) 



t;iU3 + TO2(A3 - Ai) V2V3 + ivi{X2 ~ X3) Xs'^ 



V2 



Vl 



(73) 



From the previous arguments we know that R has a real non-negative spectrum. We denote the eigenvalues of R by 
p\, P2^ /i§, and assume pi> p2> P3> ^ (again leaving apart possible degeneracies) and define 



Xi=pi-p2-P3, X2^P2-Pl-P3, X2 ^ P3 - fJ-l ~ P2, Xi = Pi + P2 + P3- (74) 
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Hence X2 = — (/ii — fJ'2) — fJ-a < 0, = — (mi — fJ-s) — tJ-2 < 0, and x^ > 0. It follows thus that the polynomial 

P{x) = {x — xi){x — X2){x — X3){x — X4) (75) 

has at least two negative roots X2 and 2:3 and at least one positive root 0:4 . The sign of the remaining one xi depends 
on Xi and Vi. 

Straightforward calculations give 

P{x) =x'^ ~ (4Tri?)a;2 _ 8(det Rf^x + 2 Tt{R^) - (Tri?)^ (76) 

If 2 Tr(i?^) — (Tri?)^ < then the consecutive signs of the coefhcients of P{x) read H , hence from the Descartes 

rule of signs it has at most one positive root and consequently xi — fii ~ IJ-2 ^ 1^3 < 0. On the other hand for 

2 Tr(i?^) — (Tri?)^ > the signs of the coefficients of P{~x) arc H hence P{—x) has at most two positive roots, 

i.e. P{x) has at most two negative roots and, consequently, xi is positive. (The validity of this assertion can be easily 
established by considering the graph of P{x) under the conditions that it crosses the x-axis in at least two negative 
and at least one positive point). Summarizing: the sign of fJ-i — H2 ^ Ms coincides with the sign of 2 Tr(i?^) — (Tri?)^, 
or in other words, p is 'classical' if and only if 2 Tr(i?^) — (Tri?)^ < 0. (The equality in the last formula appears when 
we relax the non-degeneracy conditions). 

It can be proved that ([7T|) is equivalent to the criterion of [ij for the classicality of spin-1 coherent states. The 
latter is formulated as follows. We rewrite p in the form 

P=\{I-W + u-L), (77) 

where u is a real vector and L = (Li, L2, L3) with Li given by ((67)) . The matrix Z of Giraud et.al. is then defined as 
Zij = Wij — UiUj. The state is classical iff Z is positive definite. The proof of the equivalence of both criteria is given 
in the Appendix. 



VI. 'CLASSICAL' STATES AS GENERALIZED COHERENT STATES 



As already stated in the introduction similarities of all above outlined constructions have their common origin in 
the fact that 'classical' states are special orbits of the Lie group of symmetries characteristic for a particular problem. 
To make this statement more precise let us recall a few facts from the theory of group representations needed for 
the definition of generalized coherent states [25|] . Let K he a. compact semisimple Lie group, t its Lie algebra and 
Q = l + il- the complexification of I. As a convenient example we can take K = SU{N), i = sun, and g = sIn{C). In 
addition we can imagine all of them as sets of complex matrices in the defining representation in the complex space 
C^, then K = SU{N) is the set of unitary N x N matrices with determinant 1, t = suat — the set of traceless 
antihermitian N x N matrices, and 51n{C) — the set of complex, traceless N x N matrices. 

The Lie algebra g can be decomposed into the direct sum 

g = n_eten+, (78) 

where t is the Cartan subalgebra of g (its maximal commutative subalgebra dimension of which is called the rank of 
g), whereas n± are particular nilpotent subalgebras of g. In our SU{N) example realized in we can choose t as 
the set of traceless diagonal matrices and n+ and n_ as, respectively, upper- and lower-triangular matrices. 
As linear spaces n-(- and n_ are direct sums of the root spaces 

n+ = 0g<„ n_=0g_„. (79) 

a a 

Each go, is one-dimensional and is uniquely determined by the commutation relations of its elements with the elements 
in t 

[H, X] = a{H)X, H et, X G g„, (80) 

where a(-) is an appropriate linear form on t. If we choose some basis {Hi}, {i — 1, . . . ,r = dimt) in t, we can calculate 
(|80| for H = Hi obtaining vectors (roots) a with the components — a{Hi). They span the Euclidean space of 
dimension equal to the rank of g, in which we can choose a basis consisting of positive simple roots — each root a 
is a linear combination of them with only non-negative (positive roots) or non-positive (negative roots) coefficients. 
The algebra g is uniquely determined by its positive roots (or, equivalently, by its positive simple roots). The positive 
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(negative) root vectors are, by definition, the elements X±a, of Q±a, fulfilling, according to ((80)l . [Hi,X±a] = ±aiX±a. 
There is one-to-one correspondence between the negative and positive root vectors. 

On a semisimple Lie algebra g there exists a nondegenerate bilinear form (the Killing form) defined as 

{X, Y) = Tr(adX • adF), (81) 

where adX is the linear transformation of g into itself given by a,dX{Y) = [X, Y]. We can use the Killing form to fix 
the normahzation of root vectors and the elements Hi by choosing {Xa, X-a) ~ 1 and {Hi, Hj) — Sij. 
From the root vectors and the elements Hi we can construct the second order Casimir operator 

T 

C2 ^(X„X_„ + X_„X„) + Hf, (82) 

a>0 1=1 

which, as straightforward calculations show, commutes with every element of g. With the help of C'2 one can conve- 
niently characterize the set of coherent states. 

The group K as well as the algebras t and g can be irreducibly represented not only in the defining space C^, as 
we did with SU{N), but also in spaces of other dimensions. If we choose as the representation space H = C*^, we 
can again think of representatives of K, g, etc. as sets of matrices, so we will use the same letter tt to denote the 
homomorphism between the sets of matrices in and C*^ defining the considered irreducible representation. Thus 
we denote by Tr(H) the matrix representing in C*^ the element H of, say, g, or by Tr{U) a representative oi U E K 
etc. 

For each irreducible representation of a semisimple g there exists a vector \vmax) m Ti which is a common eigenvector 
of all n{Hi) annihilated by 'k{X) for X E n+, i.e. 

'^iHi)\v,nax) = Aj|Wmax), 7r(X)|Wmax) = 0, for X G n+ . (83) 

The vector \vmax) is called the highest weight vector of the irreducible representation. An irreducible representation 
of g, as well as K, is uniquely determined by the vector A = (Ai, A^). 

Generalized coherent states for the group K, or more precisely, for its irreducible representation in 7i, are obtained 
by applying to the highest weight vector \vmax) all possible representatives 7r(/7), U £ K (s^. In all cases relevant for 
the present considerations, we can assume that the representation in question is unitary, so the action of the group 
via its representation does not change the lengths of vectors in Ti., and we can assume \vmax) to be normalized to the 
unit length. Since in quantum mechanics we identify states differing by a phase factor, it is more convenient to look 
at the coherent states as an orbit of K in the projective space FTC, ie. the space of the equivalence classes of vectors 
differing by a multiplicative constant. Let us denote the equivalence class of some |w) S 7i by [{w}] 

[\w)] = {\w') e n : \w') = c\w), c e C}. (84) 

Since the representation tt acts in a natural way in P7i, 

AU)[\w)] = k(C/)k)], (85) 

we may identify the manifold of coherent states Cs with the orbit of K in P7Y through the equivalence class (the 
complex line) of the highest weight vector [Ifmaa;)] 

Cs = {Ti{U)[\v„,ax)]-U eK). (86) 

Equivalently, we may consider Cs as a manifold K/K^, where is the isotropy group of \vmax), 

K, = {UeK: 7T{U)\v^ax) = e'^lv^ax), for some a e R}. (87) 

The isotropy group of the maximum weight vector [Iwmaa;)] depends not only on the group K, but also on the 
representation tt. 

The imaginary part of the scalar product in a complex Hilbert space Ti. defines a linear symplectic structure in 
it (or more precisely in the real space obtained from Ti. by treating real and imaginary parts of the components of 
vectors in Ti as independent real components of vectors in a real vector space) . This symplectic structure is inherited 
by FTi. The manifold ((86|l treated as a submanifold in PTi happens to be a symplectic one, ie. can be used as a 
classical phase space. In addition it is a complex manifold (in fact a unique one among the orbits of K in PTi), which 
makes it a Kahler manifold [26j . 
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The set of coherent states Cs, ie. the orbit through the highest weight vector can be uniquely characterized by 
certain bihnear condition Let us denote by 6 the vector which is the half of the sum of all positive roots a. A 
state \u) is coherent if and only if 

^ = (2A + 2d, 2A) ® , (88) 
where (, ) is the scalar product in the r-dimensional Euclidean space of vectors a and d, and 

r 

L = 7r(C2) (E)I + I(E) 7r(C2) + ^ (tt{X^) ® 7r(X_„) + 7r(X_„) ® 7r(X„)) + 2 ^ 7r(iJ,) » 7r(i7,)- (89) 

Q>0 1=1 

We used the above formula, adapted to the case of SU{2), to write the bilinear characterization of spin-coherent states 
(|65|). however the other bilinear characteristics of 'classical' states (cf. (|34|) . (ISTl and ((6T|) ) where not straightforwardly 
derived from (j89p . The connection between two different bilinear characterizations will be explained elsewhere. 

All classes of pure states considered in the preceding sections can be treated as particular instances of the general 
coherent states. Thus 

1. spin j coherent states are generalized coherent states for the SU{2) group irreducibly represented in C^-'"'""'^, ie. 

K = SU{2), H = C^^+\ (90) 

2. for the fermionic separable states the appropriate identifications read 

K = SU{N), H^A{C^(E)C^), (91) 
where A denotes the antisymmetric part of the tensor product, 

3. for the bosonic separable states we have 

K = SU{N), H = S{C^ (E)C^), (92) 
where S denotes the symmetric part of the tensor product, 

4. for the separable states of two distinguishable particles, 

K = SU{N) X SU{M), H^C^ (E) C^'^ (93) 

VII. LOW DIMENSIONAL 'CLASSICAL' STATES AND ORTHOGONAL SYMMETRIES 

Low-dimensional states considered above, apart from the exhibited properties of being effectively characterizable by 
concurrences based on antilinear operators, have some additional nice representations which we would like to outline 
shortly in the present section. 

Let us start from the following observation. A point in complex projective space of an arbitrary finite dimension 
(i.e. a pure state of a quantum system with an iV-dimensional Hilbert space Ti.) can be always represented in the form 

IV') = cos^x + isinfly, (94) 

where x and y are real, unit, orthogonal vectors 

x2 = y2 = 1, X • y - 0. (95) 

Indeed, an arbitrary vector Cz Ti can be decomposed into its real and imaginary parts, = x' + iy' . The 
normalization condition 1 = = x'^ + y'^ leads to x' — cos^x, y' = sin6'y with real, unit vectors x and y. 

Since we work in the projective space FTC, we still have at our disposal a phase factor, \ip) ~ e^'^jV'), which, as simple 
calculation shows, is enough to enforce orthogonality of x and y. Obviously, by alternating the global signs of x 
and/or y, we may restrict 6 to the interval [0, j]. 

The representation (j94p . (|95[) is, in general, not particularly useful or interesting, since it is not invariant under 
general unitary transformations of 7i, which are allowed by quantum mechanics, but in general do not leave x and 
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y real. If, however, for some reasons the relevant transformations made on the system in question belong to the 
orthogonal group SO{N), the situation changes. 

For all low dimensional cases considered in the previous sections this is the case. For two qubits, although the 
Hilbcrt space is four dimensional, with the 'natural' symmetry group S'[/(4), the relevant subgroup of the local 
unitary transformations which do not change the entanglement is SU{2) x SU{2), homomorphic to the orthogonal 
group SO{4). For two spin-3/2 fcrmions the Hilbert space is six dimensional but the relevant subgroup of SU{6) of 
transformations which respect the antisymmetry of states is SU (4) — the full symmetry group of the four-dimensional 
single-particle space (the antisymmetry is not destroyed if both particles undergo the same operation). Here again 
we have a homomorphism with an orthogonal group: S'[/(4) ^ 5*0(6). Finally, in the case of two bosons in three- 
dimensional space (two-dimensional single particle space) and coherent spin-1 states the relevant subgroup of SU{3) 
is SU{2) — the full symmetry group in the single-particle space in the case of bosons and the group defining the 
structure of coherent states, and the relevant homomorphism is the familiar one SU(3) ~ 5*0(3). 

Throughout we consider the case when a group of transformations leaves the set of 'classical' states invariant (or 
more generally, when the generalized concurrence is constant on the orbits of the group). A recurrent theme is that 
the group is, on the one hand a unitary subgroup SU{K) of the special unitary group SU{N) of the underlying 
Hilbert space, and on the other hand it is homomorphic to the special orthogonal group SO{N). The examples are 
exhausted by the following table: 

SU{6) D 5C/(4) - 50(6) (spin-3/2 fermions), 

SU{4) D SU{2) X SU{2) - 50(4) (qubits), (96) 

SU{3) D SU{2) - 50(3) (2-dim bosons, spin-1 coherent states). 

There are no more homomorphisms of the above type, as it is clear from the inspection of the existing isomorphisms 
between unitary and orthogonal Lie algebras (c.f. [2g|, p. 519-520). This observation underscores again the exceptional 
features of low-dimensional systems. 

To take direct advantage of the representation ((94)) one should determine the basis in which the relevant group is 
represented by real matrices. A moment of reflection suffices to establish that such a basis consists of vectors fulfilling 
I'^fe) = ©I'^fe)- Indeed, in this basis a unitary transformation i— > U\ip) leaving (V'l jSIV's) invariant, fulfills U^U = /, 
i.e. U is orthogonal (and unitary, hence real). Such bases are dubbed 'magic' — for two qubits the magic basis consists 
of four orthogonal Bell states [l3|, in the fermionic and bosonic cases they are given in [l^ . 

The angle 9 and vectors x and y for 'classical' states must fulfil {iIj\T\iIj*) = (cosfe— i sin0y)-T-(cos 9x—i sin By) ~ 0. 
For spin-1 coherent states this gives particularly simple condition cos^ 9 — sin^ 9 = 0, i.e. 9 = n/A. 



VIII. ACKNOWLEDGMENTS 



The support by SFB/TR12 Symmetries and Universality in Mesoscopic Systems program of the Deutsche 
Forschungsgemeischaft and Polish MNiSW grant no. DFG-SFB/38/2007 is gratefully acknowledged. IB was sup- 
ported by the Swedish Research Council, VR. 



IX. APPENDIX 



We will prove that positivity of the matrix Z defined by Giraud et al. is equivalent to 2Tr(i?^) — (Tri?)^ < 0. For 
p given by ([72]) we get 



Z 



-Al + A2 + A3 - 4?;i^ -4:V2Vi -ivsvi 

-4:V2Vi Al - A2 + A3 - 4U2^ -4:V2V3 

-4u3Ui -4w2«3 Al + A2 - A3 - 4U32 



(97) 



After a short calculation we establish that the coefficients of characteristic polynomial Pz{x) = -f a2X^ -f aix + oq 
of Z read 



02 = 4{vl+vl + vl)-l, 

fli = Bdetp + 4(AiA2 -I- A1A3 -I- A2A3) - 8A1A2A3 - 1, 
ao = 2Tr(i?2) - {TiRf. 



(98) 
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We will prove the following 

Lemma 2 Let: ri > 0, Si > 0, i = 1, 2, 3 and 



7-1 + ra + = 1 (99) 
rirarg > risi + r2S2 + rgsa. (100) 



Then 



Sl+S2+S3<^. (101) 

Proof 

Let, eg., 1 > ri > r2 > ra > 0. Then 

rir2r3 > r3(si + S2 + S3) + (''i - ''3)si + (^2 - r3)s2 > ^3(51 + S2 + S3). (102) 

Hence 

''l''2 > Si + S2 + S3. (103) 

From the arithmetic-geometric mean inequality 

V^<'-^<1, (104) 

hence 

1 , , 

- > rir2 > si + S2 + S3. (105) 

□ 

We substitute in the above lemma = and Si — vf. Then (|102p and (|103p reduce, respectively, to the true 
statements Tip = 1 and detp > 0. As a consequence we get 02 < 0. We have thus the following possibilities 

1. ao < 0, ai > 

In this case the signs of the coefficients of Pz{x) are H 1 — , whereas those of Pz{—x) read , which 

means that all three roots of Pz{x) are positive, hence Z is positively definite. 

2. ao > 0, ai > 

The signs of the coefficients of Pz{x) read H (-+, i.e. Pz{x) has at most two positive roots. Consequently, at 

least one root is negative and Z is not positively definite. 

3. flQ > 0, fli < 

The signs of the coefficients of Pz{x) read -i h which again gives at most two positive roots, hence Z is not 

positively definite. 

The remaining case ao < 0, ai < is excluded. Indeed, we will prove ao < ai > 0. To this end assume as 
previously Ai > A2 > A3 > and define: 

(71 = A2 + A3 - Ai, g2 = Ai + A3-A2, g3 = Ai + A2-A3. (106) 

The implication we want to prove reduces now to 

91(7293 > 4qi(72W3 + 4qi(73W2 + 4(72q3Wi (107) 

(7192+9193 + 9293 > 4((7i + 52)^1 + 4((7i + (?3)w2 + 4(92 + 93)w?- (108) 

We have 91 + 92 + 93 = Ai + A2 + A3 = 1 and (72 > 0, ^3 > 0, whereas the sign of qi can be arbitrary. 
If gi > we have from ([TU7| 

919293 > 49i92w| + Miq^vj + 4(7293^? > 49193^2 + 49293^?, (109) 
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hence 



9192 > 4(7iW2 + 4g2Wi- 



(110) 



Analogously 



(111) 



and 



92(73 > 493 ^^2 + 492^^1- 



(112) 



Adding (|110p - (|112|) we obtain the desired result. 
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